Drastically suppressing the error of ballistic readout of qubits 
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The thermal jitter of transmission of magnetic flux quanta in long Josephson junctions is stud- 
ied. While for large-to-critical damping and small values of bias current the physically obvious 
dependence of the jitter versus length a ~ y/L is confirmed, for small damping starting from the 
experimentally relevant a = 0.03 and below strong deviation from a ~ y/L is observed, up to nearly 
complete independence of the jitter versus length, which is exciting from fundamental point of view, 
but also intriguing from the point of view of possible applications. 
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In the recent years superconducting circuits have at- 
tracted a considerable interest as promising devices for 
quantum computations jll-Ql ■ The advantage of supercon- 
ducting circuits in comparison with other types of qubits 
is the possibility to combine in one chip both the qubits 
and the readout electronics. Mostly, dc SQUIDs are used 
for the readout from superconducting qubits. Recently 
it has been suggested to use for the readout the well- 
elaborated Rapid Single Flux Quantum devices (RSFQ), 
redesigned in a special way. The main idea of this type 
of readout is the ballistic propagation of magnetic flux 
quanta along two separated Josephson transmission lines 
(JTL) Q-^j consisting of either one long or series of 
short junctions having weak damping, see Fig. [TJ If one 
JTL is inductively coupled to the qubit, the magnetic 
field of qubit will change the speed of the fluxon (soli- 
ton), traveling from one JTL end to another one, and 
thus will change the fluxon transmission time. This dif- 
ference of fluxon transmission times in two JTLs can ef- 
fectively be measured using the existing RSFQ circuitry 
(SFQ receiver in Fig. 1), while the proper fluxons (SFQ 
pulses) can be produced in SFQ driver. The use of un- 
derdamped JTLs together with the standard RSFQ cir- 
cuitry in close vicinity to a qubit is very promising, since 
allows to eliminate expensive external pulse generators 
and readout equipment, to nearly eliminate parasitic ca- 
pacitances and inductances, altering the readout pulses, 
and to finally improve the readout fidelity of the qubits. 
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FIG. f : The illustration of the ballistic readout from the qubit 
using two Josephson transmission lines, where top one is in- 
ductively coupled to the qubit. 

It is obvious that increasing the JTL length, the pre- 



cision of time difference measurement can be improved. 
However, different types of fluctuations restrict the read- 
out precision since they lead to the jitter of traveling 
fluxons. The importance to study jitter in Josephson 
junctions was first understood in Q, see also Later 
jitter in short junctions has been studied both analyti- 
cally and numerically [ll|- 14 1. However, as it has been 

llij . for se- 
N, where 
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understood and experimentally verified 
ries of junctions the jitter increases as a 
N is the number of junctions. Therefore, for JTLs the 
jitter must increase as cr ^ \/L, where L is the JTL 
length, which also follows from simple theory presented 
in [7]. Since the fiuxon traveling time r is proportional 
to the JTL length (r ^ L), one must find a compromise 
length, leading a low readout error a/r. However, to 
our knowledge, except simple reasoning and the theory 
of [7], the jitter in underdamped JTLs and long Joseph- 
son junctions have not been studied neither theoretically 
nor experimentally. The experimental results devoted 
to RSFQ can not be directly applied since they are ob- 
tained for critically damped junctions rather than for un- 
derdamped. The present paper, therefore, is devoted to 
detailed theoretical study of the jitter versus damping 
and the length in an underdamped long Josephson junc- 
tion. While for large-to-critical damping and small values 
of bias current we fully confirm a ^ \fL dependence of 
the jitter versus length [7], for small damping starting 
from the experimentally relevant a = 0.03 and below, 
we observe strong deviation from a ^ \/L, up to nearly 
complete independence of the jitter versus length, which 
is exciting from fundamental point of view, but also in- 
triguing from the point of view of possible applications. 
In particular, this allows to drastically suppress the read- 
out error by increasing the length of underdamped JTLs. 

In the course of the present paper let us consider the 
long Josephson tunnel junction (JTJ) as a limiting case 
of JTL with very little spacing between short junctions, 
much smaller than the Josephson penetration depth Aj. 
All calculations are performed in the frame of the sine- 
Gordon equation, giving a good qualitative description 



of basic properties of long JTJ: 

<j>tt + OKpt - ipxx^i- sin((/)) + if{x, t), 



,{t) = 



dP{t) 



(1) 



where indices t and x denote temporal and spatial deriva- 
tives. Space and time are normalized to the Josephson 
penetration length Aj and to the inverse plasma fre- 
quency a;~^, respectively, a — ujp/ujc is the damping 
parameter, ujp = •\/2e/c/?iC, Wc = 2eIcRN Ic is the 
critical current, C is the JTJ capacitance, i?jv is the nor- 
mal state resistance, i is the dc overlap bias current den- 
sity, normalized to the critical current density Jc, and 
if{x, t) is the fluctuational current density. If the critical 
current density is fixed and the fluctuations are treated 
as white Gaussian noise with zero mean, its correlation 
function is: {if{x,t)if{x',t')) — 2a'-fS{x — x')5{t ~ t'), 
where 7 = It/{ Jc^j) is the dimensionless noise intensity 
17| . It = 2ekT/h is the thermal current, e is the elec- 
tron charge, h is the Planck constant, k is the Boltzmann 
constant and T is the temperature. 

The boundary conditions that describe coupling to the 
environment, have the form (isj : 

0(0, t)x + rLCL^O, t)xt - CL0(O, t)tt - r, (2) 
t)x + rRCR(f>{L, t)xt + CR(f>{L, t)tt = r. (3) 

Here F = is the normalized magnetic field, and L is the 
dimensionless length of JTJ. The terms with the dimen- 
sionless capacitances and resistances, clu and r^^/j, are 
the RC-load of a JTJ placed at the left (input) and at 
the right (output) ends, respectively. 

For simplicity we assume that the bias current density 
is uniformly distributed along the space i{x) — const. As 
initial condition a kink (j){x,t) — 4 arctan [exp(a; — xq)] 
inside the junction is taken (for xq = 5). Temporal 
and spatial intervals are chosen Ai = Ax — 0.01, and 
it has been verified that further decrease of the steps 
did not change the results. Two values of i?C-load were 
tested: mismatched case where simple boundary condi- 
tions d(j){x = 0)/dx = d4>(x = L)/dx — were used, and 
perfectly matched case = rn = 1, cl — cr = 100. 
Since the results are nearly the same, the curves for the 
perfectly matched case (to suppress reflected waves from 
the junction ends, as required in experiments) are shown. 

Two different definitions of the fluxon traveling time r 
and the jitter a are used. One is the usual mean first pas- 
sage time of the boundary, where random time at which 
the soliton hits the right junction end is computed and 
after the mean value and the standard deviation are cal- 
culated by averaging over 1000 realizations. Another def- 
inition is based on the notion of the integral relaxation 
time [1^ , U 1 , 12 h if the probability to find the soliton in- 
side the junction is computed P{t) = W{x,t)dx, the 
mean traveling time and the standard deviation (jitter) 
are 



dt[P {00) - P{0)]' 



Since we have investigated our task in the limit of small 
noise, both definitions give the same results within cal- 
culation precision as it must [l9j . 
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tw{t)dt, a = ^{t^) - {tf, (4) 
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FIG. 2: Jitter versus JTL length for various values of damping 
a and bias current i for 7 — 0.001. Symbols - simulations, 
dashed curves - theory (|4)). (|10p . 

In Fig. [3] the plots of the jitter versus JTL length are 
shown for various values of damping a and bias current i, 
and for noise intensity 7 = 0.001. One can see that while 
for moderate values of damping a = 0.5 and a = 0.1 the 
curves are perfectly fitted by the dependence a ~ \/L, 
starting from a = 0.03 the jitter is proportional to cr ~ 
v^Z and for decreasing values of damping and increasing 
bias current the jitter tends to a constant, which is the 
main and quite surprising result of the present paper. 

In order to understand the origin of the observed result 
let us consider the evolution of soliton shape. From Fig. 
|3]one can see the snapshots of the fluxon accelerating and 
contracting during its motion along JTL. Analyzing the 
soliton dynamics for different values of damping and bias 
current it has been noted that if the damping is moderate, 
and the bias current is small, then the soliton's shape 
weakly changes with time. In this case the dependence 
of the jitter is well described by u ~ \/L. However, if 
the damping is small, see, e.g.. Fig. [3] for i = 0.1 and 
a — 0.03, the Lorentz contraction is observed, the soliton 
width seriously decrease near the output boundary in 
comparison with the input one. Considering this as the 
main explanation that can be given for violation of tr ~ 
\/L dependence of the jitter, we have tried to obtain more 
advanced theory, since the formula of ,7] does not take 
into account the nonstationarity of the fluxon dynamics. 

Starting from Eq. ([T|), one can write the equation for 
the momentum p{t) = — | (j)x4>tdx of the center of 
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FIG. 3: Snapshots of the sohton, accelerating and contracting 
during the motion along JTL, i — 0.1 and a = 0.03. 



the soliton [20| with account of effect of noise 



dp 1 



(5) 



where p = v/ \/l ~ v^, v{t) is the fluxon velocity, and the 
noise intensity (C(OC(^')> = ^JT-v^ ~ "^^^^ ^ 
same as in Eq. ([T|). Directly from Eq. ([5]), neglecting 
the effect of noise, one can derive the expression for the 
fluxon velocity v{t) for u(< = 0) = 



v{t) 



1 



(6) 



From Eq. ([5]) one can derive the following equation for 
the location of the center of the fluxon X{t) 



X{t) = -aX{t){l 



') + {\7r^+m)va-v'r, (7) 



where X{t) and X{t) stand for temporal derivatives. 
This equation is too complex to be solved analytically, 
and we will use the following approximations. First, let 
us neglect random component of velocity (which can be 
done in the limit of small noise), so v{t) is described by 
Eq. (|6]). Second, let us also neglect by the term (1 — w^) 
in front of the first temporal derivative X{t), so Eq. ([7]) 
transforms into the equation of massive Brownian parti- 
cle, but with the noise intensity, depending on time. This 
Brownian diffusion is described by the Gaussian proba- 
bility density, since the noise ^{t) is Gaussian. Therefore, 
one can use the approach described in the chapters 1 and 
3 of the book [22], and can get the following expression 
for the variance D{t) of the process X{t), taking into 
account the nonstationarity of the noise intensity 
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2e-"^-He-2"2/] (l_„(y))5/2^y^ (8) 



where is given by Eq. ([5]). Analogically, the mean 
■m{t) of X(t) can be derived (with X{t = 0) = Xq) 



m{t) = Xo+ v{y)dy 



(9) 



Substituting the mean and the variance into the Gaussian 
probability distribution one can obtain the probability to 
find the soliton inside the junction 



P(i) = l-erfc {L-m{t))/y/2D{t) /2. 



(10) 



Substituting P{t) into the definitions (|3]), the mean trav- 
eling time r and the standard deviation a can be com- 
puted. The comparison between computer simulation re- 
sults and the theory is presented in Fig. I2|4l l6l One can 
see good agreement for larger damping, while for smaller 
damping the theory slightly underestimates the numeri- 
cal results, presumably due to our simplifications in Eq. 
([7]). Nevertheless, the agreement is good enough and 
correctly describes the qualitative behavior, confirming 
violation of \/L dependence of the jitter. 
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FIG. 4: Jitter versus JTL length for different values of bias 
current i and a — 0.01, 7 — 0.001. Symbols - simulations, 
dashed curves - theory (|4|l. (|10|l . The dependence changes 
from \/L, (top curve, circles) to \^ (diamonds). 

Let us perform calculations for fixed small value of 
damping a — 0.01 and for various bias current values, 
see Fig. |4l It is seen that for small values of bias cur- 
rent the dependence a ^ \fL is maintained, while for 
larger values starting from i = 0.01 the deviation is ob- 
served. It should be noted, however, that one must think 
about thermal budget of both qubits and its readout elec- 
tronics, trying to decrease the total heat, since usually 
cryostats have rather low thermal power at low temper- 
atures. Nevertheless, the bias current values of order 
i — 0.01, at which the deviation from the dependence 
a ~ \fL is observed, seem to be rather low to maintain 
low heat emission. 
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cr readout error a jr. 

— r- • ^ In conclusion, the thermal jitter of transmission of 

magnetic flux quanta in long Josephson junctions is stud- 

* led. While for large-to-critical damping and small values 

of bias current the known dependence of the jitter ver- 
—<>- * igng|;j2 (J ^ yfij ig confirmed, for smaller damping the 

strong deviation from a ~ y/L is observed, up to nearly 

^ x_ -X _x-x complete independence of the jitter versus length, which 

^ ^ ♦ 1=0 001 y=oooi exciting from fundamental point of view, but also in- 

1=0.002, y=o.ooi triguing from the point of view of possible applications. 
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FIG. 5: Jitter versus JTL length for very low damping a = 
0.001. Symbols - simulations, dashed curves - theory (|4)l. (|10p . 
All curves are well fitted by \fL dependence. 

It is known that with decreasing the temperature of 
JTL from 4K to 50 mK the damping decreases by one 
order of magnitude . This allows to operate with much 
smaller bias currents to decrease heat emission, see Fig. 
[5j The curves, calculated for a = 0.001, demonstrate 
good fitting to the dependence tr ~ even for very low 
bias current i = 0.001. 
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FIG. 6: The readout error cr/r versus JTL length. Symbols - 
simulations, dashed curves - theory (|4|). (|10p . The dependence 
changes from 1/\/L (top) to l/L (bottom). 

Finally, let us consider the ratio of jitter to the travel- 
ing time (t/t, which determines the error of measurement 
(readout error) due to thermal fluctuations. In Fig. |6]the 
plots of (t/t are presented versus JTL length, for differ- 
ent parameters a and i. One can see that contrary to the 
conclusions of Ref. ^] larger bias current leads to smaller 
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